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O Abstract 

.The present work is a reply to the paper [ 1 ] . It is proven that the argumentation of Ref . HI is inconsistent. The variational functional 



^ for the polaron ground state energy considered in Ref. [1] contains an incomplete recoil energy. Since the variational functional 
of Ref. U] is incomplete, it is not proven to provide a variational upper bound for the polaron ground-state energy. The same 
conclusion follows also for the bipolaron ground-state energy. 
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Polarons and bipolarons are invoked in the study of polar ma- 
terials, including high-r t . superconductors jp, La |4J]. Rigorous 
variational methods (see, e. g., Ref. LTD) are important in 
this field, i. a. because in the bipolaron mechanism of super- 
conductivity, the parameters of the superconducting state and 
the critical temperature strongly depend on the bipolaron bind- 
ing energy. 

The work (jl is a reply to our comments JH] on the variational 
approach aimed at in Refs. $M\ 1 111 . In Ref. QSQ we show that 
the strong-coupling expression for the bipolaron ground state 
energy calculated in Refs. JslQli is not justified as a variational 
upper bound. 

It is suggested in Ref. [1] that a properly chosen cutoff for 
the phonon momenta leads to correct variational polaron and 
bipolaron ground-state energies in the strong-coupling limit. 
However, this conclusion is not valid, because the recoil en- 
ergy treated in Refs. J9, 1(1 HI is incomplete, as we wrote in 
Ref. H. 

The complete polaron recoil energy within the approach of 
Ref. [Ill] was found by Porsch and Roseler [ 12]. They showed 
that, when imposing a cutoff for the phonon momentum, the 
polaron recoil energy Er consists of two parts: 
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where EjP is the recoil energy determined in Ref. | IT], and the 
term dE { R PR) is given by Eq. (43) of Ref. JH: 
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where qo is the cutoff value for the phonon momentum, a> q = 
a>o + with a>o the LO-phonon frequency, and lo. q } are the fre- 



quency eigenvalues resulting from the Bogoliubov-like canoni- 
cal transformation for the phonon operators (performed in Refs. 

EMII). 

It is stated in the reply |lt] that the reasoning of Ref. lUt] 
is "based on the erroneous approach ... to the strong cou- 
pling limit when the cutoff parameter is introduced in the the- 
ory." However, the argumentation of Ref. QJJ] is related only to 
the term E ( p, ignoring the Porsch — Roseler term SE R \ In 
the present work we treat the contribution to the recoil energy 
dE R PR) missed in Refs. M, HI. 

The expression obtained in RefJll for Q. qo readsS 
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with the function 
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Here, f (q) are variational functions. In Ref. 
chosen as 
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ll ill , they are 
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with the variational parameter a and the amplitudes of the 
electron-phonon interaction V q . 

In Fig. 1, we plot the complete recoil energy Er and the 
contributions Ei T , dE^^ as a function of a for qo — 8 and 
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'There is a misprint in Eq. (42) of Ref. 1 12] corresponding to Eq. (5) of the 
present work: the factor 3 in the denominator is missing. 
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a — 4 (measured in units of J^jf-)- The arrow indicates the 
value of the coupling constant 



at which the steep maximum of the integrand in E R 
tioned in Ref. II ill t crosses the cutoff boundary. 
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For sufficiently small a, the Tulub's recoil energy E R dom- 



inates, and 6E R PR> is negligibly small. When a increases (keep- 

( T) 

ing other parameters constant), E R tends to a finite value, 
while 6E R monotonically increases. 
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Figure 1: The recoil energy En (solid black curve), the contributions E R 
(dashed red curve) and SE^ R> (dotted green curve) as a function of a for qo = 8 
and a = 4. The dot-dashed blue curve is the recoil energy without cutoff |8]. 
The arrow indicates the value a c at which the peak discussed in Ref. II ill passes 
the cutoff boundary. 



The analytic asymptotics for the Tulub and Porsch — Roseler 
terms of the polaron recoil energy in the strong-coupling limit 
gives us the results 



(as in Ref. [11]), and 
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which remarkably coincides with the asymptotic strong- 
coupling expression from Ref. [8] for the recoil energy without 
a cutoff. For an increasing a, the optimal value of the varia- 
tional parameter a increases. Thus the Porsch — Roseler con- 
tribution dominates in the recoil energy in the strong-coupling 
regime, while the Tulub contribution E R T> constitutes only a 
residual part of the recoil energy. Moreover, the analytic for- 
mula ([8]l confirms the results of our previous treatment [8]. The 
same conclusion is valid for the bipolaron ground-state energy, 
because the recoil contributions in the polaron and bipolaron 
problems are structurally similar to each other. 



In addition, there are logical inconsistencies in the argumen- 
tation of Ref. lUt] (which are of a secondary importance with 
respect to the question discussed above). It is stated in Ref. 
(jl that "They ... used the asymptotics q(l/A) = 2y/3A as 
the basis for their calculations of the polaron energy." However, 
the calculations in [8] are performed on the basis of the com- 
plete expression for the function q (1/A) given by Eqs. (2.11) to 
(2.12) of the work by Tulub ifTHl rather than its strong-coupling 
asymptotics. Also it is written in QJJ] that Tulub's choice of the 
variational functions / (q) given by Eq. © is the best, while 
the choice of / (q) in Ref. [8J] is the worst. However, the varia- 
tional functions used in Ref. li are the same as those in Refs. 

In conclusion, when accounting for a phonon cutoff, the 
Porsch — Roseler contribution to the polaron and bipolaron 
recoil energy dominates in the recoil energy in the strong- 
coupling regime calculated. This contribution is missed in Refs. 
lIU |Tol [Till . As a result, the variational functionals for the 
olaron and bipolaron ground-state energies derived in Refs. 
[Tol [Till are incomplete. Consequently, these variational 
functionals are not rigorously proven upper bounds. 
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